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, I briefly review the theory and phenomenology of bottom production at colliders. 

1 When all theoretical uncertainties are taken into proper account, and modern 

' measurements are considered, no significant discrepancy is seen between data and 

' QCD predictions 

(N : 

I The physics of bottom quarks is one of the best-studied topics in particle physics. 

^ I Experimentally, this is due to the abundance with which h quarks are produced at 

i colliders. Theoretically, the reasons of interest are many. The characteristic that 

' sets the b apart is its heaviness, 

cn ; 

^ i TOfc > Aqcb (1) 

\l ' which entails peculiar properties. If one is interested in the phenomenology of the 

. decays of the 6-flavoured hadrons, eq. (^) suggests to treat the b as infinitely heavy 

' in comparison with its companion light quark(s) in a bound state, paving the way to 

HQET and its symmetry properties. On the other hand, if one aims at studying the 
^ hard production mechanism, eq. 1^ implies the possibility of computing the open-fe 

O . cross section, which is free of coUinear and infrared singularities order by order 

' in perturbation theory (as opposed to, say, open-w cross section, whose final-state 

^ . collinear singularities are cancelled only upon convolution with a non-perturbative 

fragmentation function). The bottom is also the heaviest quark which hadronizes 
before decaying, allowing us to test many of the ideas of the factorization theorems 
' in a relatively clean environment. 

To be definite, let me consider the b cross section in hadronic collisions, which 
can be written as follows: 

d<^H,H.^,-biS) = y[ dx,dx2f^'''\x,)ff^\x2)da^,,-Js = x,X2S), (2) 



X 



and can be readily extended to other types of colliding particles. Here, fj;^^ are 
the parton distribution functions (PDFs) , and da^ .^g^ are the short-distance cross 
sections, the only pieces in eq. (0) that can be computed in perturbation theory. 
The LO term (of 0(q;|)) is trivial to obtain. The NLO one (of 0{al)) was the result 
of landmark calculations |ll2l3j . Not surprisingly, the NNLO term is not available 
at the moment; this may be worrisome, since at the NLO the scale dependence is 
still pretty large, and the corrections are 100% of the Born term. However, there 
are at least a couple of points that deserve more immediate attention than the lack 
of the NNLO contribution. The first is that, as in any other cross section computed 



© Institute of Experimental Physics SAS, Kosice, Slovakia 



1 



2 



Stefano Frixione 



in perturbation theory, large logs can appear which spoil the "convergence" of the 
series. In other words, 



da = V a,a; , a, = V a:f—>\o^'^~'^ Q , (3) 



i=2 fc=0 



where Q generically indicates a "large" quantity, in the sense that ag log^ Q^l- 
The second problem is that, although theoretically well defined, the open-& cross 
section is not physically observable. In order to compare theoretical predictions 
with data, two things can be done, a) Hadron- level experimental data are deconvo- 
luted, and presented in terms of parton-level "measurements" that can be directly 
compared with open-6 results. These deconvolutions are typically performed by 
means of parton shower Monte Carlos by the experimental collaborations, b) The 
open-& cross section is convoluted with a non-perturbative fragmentation function 
(NPFF) jj^^Hb . jtqj, ^Yjc single-inclusive Pt spectrum, one writes 

dam f dz^t^H.^^^^^da^^ ^^^^.^^ 



dpj, J z dpT 

jjb^Hb (describes how a b quark transforms ( "fragments" ) into a B hadron; it is not 
computable in perturbation theory but, being universal in the same sense as PDFs, 
can be fitted to data in a given type of collision (usually e~^e~) and used elsewhere. 

Common wisdom has it that neither strategy a) nor b) have been particularly 
successful, since Tevatron data (and, to some extent, SpS ones) have been shown 
to be systematically larger than NLO QCD predictions, regardless of whether they 
were presented in terms of b quarks or of B mesons. In a recent CDF paper ^ on 
B^ single-inclusive p^ spectrum, the discrepancy was quantified to be 2.9 ± 0.2 ± 
0.4. Taking these comparisons blindly, one is led to conclude that b physics is the 
problem of the SM, and offers the first glance beyond it '5'. Although this remains 
a viable possibility, it seems premature to buy it without first reassessing carefully 
all possible sources of mistakes in the past comparisons between theory and data, 
and considering the uncertainties that so-far uncalculated SM contributions can 
give. In particular, one should try to answer the following questions: 

1 ) Do large logs spoil the convergence of the series? 

Is the fragmentation/deconvolution performed appropriately? 
But before getting into this, let me point out that, although the discrepancies 
between data and NLO QCD have been quoted to be large, this is mainly due to 
the failure to incorporate properly all the uncertainties, including the theoretical 
ones which are very large. Upon doing so one realizes that, on a statistically sound 
basis, most of the data lie withing la from the default theoretical predictions, and 
very rarely the discrepancy exceeds the 2a level. The interested reader can find an 
informative discussion in ref. [H]. 

The answer to question 1) depends on the fact that the logarithms that grow 
potentially large can be divided into two classes. The first class includes those logs 
whose arguments don't depend on the observable being measured, such as 
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which are known as threshold logs (relevant when the cm. energy is not much 
larger than the quark mass), and small- a: logs (relevant when the quark mass is 
negligible wrt the cm. energy) respectively. Threshold logs are clearly not a factor 
at colliders; small- a; logs are estimated to give up to 30% effects at the Tevatron [7]. 
Thus, the overall picture would not change if these logs were properly resummed 
and matched with NLO QCD predictions. The logs belonging to the second class do 
have arguments which are directly related to the observables. For example, when 
measuring the single- inclusive Pt spectrum, the bb pr spectrum, or the bb azimuthal 
distance in the transverse plane, logs of the following arguments are generated 

Q=^, Q^Eim, (6) 

nib rrib n 

and the perturbative series may be badly behaved when Prib) 3> nib, Pt(^) — 
0, and A(p{bb) ~ n respectively. Since single-inclusive Pt spectra are routinely 
measured, large pT{b)/mb logs have a prominent role. The resummation of these 
logs, i.e. the rearrangement of the perturbative series of eq. © in the following 
form: 

J oc oo / 2 \ ^ 

---^EE^i''"ihl°g&) (7) 



dpi ^ ^ ■ ' V^^ m, 

can be achieved by using the perturbative fragmentation functions computed in 
ref. [Sj. Each term in the sum which runs over j corresponds to a given logarithmic 
accuracy (j = is LL, j = 1 is NLL, and so on), and PST stands for power 
suppressed terms, i.e. terms which vanish in the limit m/p^, 0. Typically, 
eq. O is computed up to the NLL; since PST are neglected, the quark behaves 
as if massless, which implies that predictions based on eq. Q must strictly be 
used only when pj, 3> mb- The trouble is that Pt S> mb is not a quantitative 
statement, and often resummed predictions are compared to data even for p^ ^ mb- 
It should be clear that such a comparison is void of sense, and any agreement 
between theory and measurement must be regarded as accidental. In any events, 
a matched computation (FONLL) was proposed in ref. |9j, which combines the 
virtues of the fixed order and of the resummed formulae: 

J OO / 2\*°0 / 2" 
2 , „ „3 , „2Y-^(0) ,_Pt \ , „3Y-^(1) ,_Pt 



^ - a,ai + a,ai + ^ r-" [as log ^ j + < 1^ r^' [c^s log ^ j (8) 

and can be used to get sensible predictions in the whole Pt range. Armed with 
eq. (IHJl, one can answer question 1) 'm a quantitative way; it turns out that, in the 
Pt range probed at the Tevatron, the effects are moderate. As discussed in ref. [TU). 
the NLO cross section used in ref. 0] is only about 20% lower than the FONLL 
one. 

Let me therefore consider question 2), and for simplicity discuss the case of 
fragmentation rather than that of deconvolution. The master equation is Q . The 
crucial point is that, while the l.h.s. of this equation is a measurable quantity, 
neither of the two terms on the r.h.s. is measurable. This is easy to understand 
if one considers that eq. Q is used to extract the NPFF from data; the l.h.s. is 
measured, the short-distance cross section da is computed, and eq. I@J is solved for 
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Figure 1. Comparison of CDF Run II data E2I with MCONLO I14I1SI and FONLL il predictions. 

NPFF: symbolically, NPFF=d(T(i/f,)/(i(7(6). Thus, if one uses da computed at the 
LO, the resulting NPFF will clearly differ from the one obtained by computing do- 
to, say, NLO. Notice that this doesn't contradict the universality property of the 
NPFF; this property merely states that a given NPFF stays the same regardless 
of the type of hard collisions involved. It follows that sensible predictions can be 
obtained only if the NPFF has been extracted from e+e" data using a cross section 
computed in the same approximation as that used to predict the pp cross section. By 
reconsidering carefully the fragmentation procedure adopted in ref . 0] , the authors 
of ref. pointed out that the claimed discrepancy of 2.9 ±0.2 ±0.4 turns actually 
out to be 1.7 ± 0.5 ± 0.5 (this also includes the 20% effect mentioned before, for 
replacing the NLO result with the FONLL one), i.e. data are within Icr from the 
default theoretical prediction. This finding is consistent with experimental evidence 
from DO that the inclusive rate of jets containing b quarks - a quantity largely 
insensitive to the details of the perturbative and non-perturbative fragmentation - 
agrees with NLO QCD predictions |12| . 

More evidence that things go in the right direction has been achieved thanks 
to a new CDF measurement ^Sl of single-inclusive 6-hadron spectrum in the 
central rapidity region. For the first time, these data probe the region of Pt — 0, 
where fragmentation effects play a minor role; thus, an improved agreement wrt the 
previous comparisons would support the conclusion that what we have to blame 
is our incomplete understanding of the fragmentation phase. This is in fact what 
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happens. The data show the best-ever agreement (see fig. ^ with FONLL and 
MC@NLO PIT^ . I stress that FONLL and MC@NLO are both based on the NLO 
computations of ref. , and they differ only in the resummation of logs beyond the 
leading ones, and in the treatment of the b quark hadronization and decay, which in 
MC@NLO is performed through the HERWIG ^B] cluster model. It is remarkable 
that, in spite of the differences in the treatment of the resummation of the large-p^ 
logs, FONLL and MC@NLO can be made to agree perfectly with a proper tuning 
of the fragmentation and of the clustering parameters (see ref. |17| for a discussion 
on this point and on the comparison between theory and Run II CDF data) . 

Let me now turn to b cross section measurements at HERA. A couple of years 
ago, the situation appeared to be consistent with what was observed at the Teva- 
tron, with NLO QCD predictions systematically undershooting HI and ZEUS data. 
This picture has now radically changed, as has been thoroughly documented at this 
conference QHI- AH data lie within 2a from the theoretical predictions; for the ma- 
jority of them, the agreement is in fact much better than 2a, with data basically 
sitting on top of NLO QCD predictions. It is worth noting that in the case of b 
production at HERA NLO and FONLL predictions coincide, since the trans- 
verse momenta probed are never too large. This also implies that the treatment 
of the fragmentation mechanism is not as delicate as in the case of the Tevatron 
measurements. 

In the case of HERA, the breakthrough that occurred in the last couple of 
years has been mainly due to the fact that, thanks to a much larger statistics, cross 
sections could be presented in the experimentally visible regions (rather than in the 
form of total rates), and compared to theoretical predictions obtained by applying 
the same cuts. Former experimental results always involved huge extrapolations 
from the very narrow visible regions to the whole phase space, performed with 
standard parton shower Monte Carlo's, which cannot give sensible predictions for 
small Pt's (see ref. for a discussion on this point), a region which gives the 
dominant contribution to the total rate. MC@NLO is reliable at small Pt's, but 
it is not yet available for photoproduction and DIS processes, and in any case the 
publication by the experimental collaborations of the visible cross sections is always 
the option to be preferred. 

In view of the lesson learned at the Tevatron and HERA, it is unfortunate 
that the measurements of the 'jj ^ bb + X cross sections |20I21) suffer from the 
drawbacks that prevented a fair comparison between theory and data in pp and 
ep collisions. The three measurements rely on huge extrapolations from the visible 
regions to the whole phase space, done with standard parton shower Monte Carlo's; 
the uncertainties associated with the theoretical predictions are too small; the tech- 
niques used are identical. For these reasons, I find it difficult to take at face value 
the discrepancies quoted (data are more than a factor of three larger than the de- 
fault NLO predictions), since a careful computation of all the uncertainties involved 
(for example, those relevant to the extrapolation to the whole phase space, which 
needs to be assessed by using at least two different Monte Carlos, and ideally the 
NLO computations themselves) would presumably show that data lie within less 
than 3a from theory. Clearly, this statement cannot be proved (or disproved) but 
by the experimental collaborations, which should follow the strategy set by HI and 
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ZEUS, of quoting results relevant to the visible regions. Let me conclude by men- 
tioning that I'm not aware of any beyond-the-SM mechanism, let alone higher-order 
QCD corrections, which could explain these huge discrepancies. 

In summary, thanks to improvements on both the experimental and the theo- 
retical sides, b data at colliders seem to be in fair agreement with QCD expecta- 
tions. I stress that the theoretical predictions are essentially based on the NLO 
computations of the late 80 's |1I2I3) . and that the changes in the predictions of 
single-inclusive Pt spectra are due to a better understanding of the fragmentation 
mechanism and to the use of more precisely determined PDF sets. Newly developed 
tools, such as MC@NLO, will serve to pin down discrepancies between theory and 
data in yet unexplored corners of the bb phase space. 
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